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INTRODUCTION 


The Naval Air Development Center (NAVAIRDEVCEN) has recently com- 
pleted the design, analysis, fabrication and static testing of a 
graphite-epoxy composite wing for the unmanned BQM-34E aircraft. A 
weight saving of 54 percent has been achieved while meeting all static 
strength and static and dynamic aeroelastic requirements. 


The BOM-34E is a high-altitude, supersonic, recoverable aerial tar- 
get vehicle. The wing has a 9-foot span and an aspect ratio of 2.5. 
Other dimensions are shown in Figure 1. The main panel of the produc- 
tion design consists of chem-milled stainless steel skins bonded to an 
alumimm honeycomb core. The outboard panels are made of fiberglass- 
reinforced plastic. With the exception of the outboard panels, the 
entire wing is factory-assembled by adhesive bonding and riveting. Its 
total weight is 157 pounds, about 30 percent of the structural weight 
of the aircraft. 


In the composite wing design, emphasis was placed on reducing the 
weight and the number of major subassemblies while maintaining the 
original airfoil shape and planform. The development was initiated under 
Independent Research (reference (a)) as an in-house effort whose ob- 
jectives were to exploit the improved material properties of advanced 
filamentary reinforced composites, to obtain an optimm design with 
minimm weight and to allow acquisition of service experience on a 
primary structural component of a high performance vehicle with no 
risk to human life. The final design, analysis, fabrication and testing 
of the wing were carried out entirely in-house under reference (b). 
Following the successful static testing at NAVAIRDEVCEN, flight testing 
commenced at the Naval Missile Center (NAVMISCEN), Point Mugu, Ca. 


Results of static and flutter analyses of a preliminary wing design 
have been reported previously in references (c) and (d). The static 
aeroelastic and flutter analyses and static test of the final design 
are described and their results presented in references (e) and (f). 
This report presents the procedure and results of the final stress and 
vibration analyses and includes the following: 


a, composite and core material characteristics including material 
stiffness properties and allowable stresses; 


b. descriptions of the wing structure and the finite element model 
used for the analyses; 


c. results of the stress and vibration analyses; 


d. computation of critical stresses and margins of safety for the 
skins and core. 
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DESIGN CRITERIA 


The flight conditions for the composite wing are the same as for 
the production metal wing, and are described in reference (g), in 
which the critical flight conditions and the resulting shear, torque 
and bending loads are presented. The equivalent static test loads are 
specified in reference (h). 


The exposed planform and aerodynamic profile of the composite wing 
are identical to those of the metal wing. Size and location of fuselage 
attachment bolts are the same to allow installation of the wing without 
modification of the fuselage. In addition, the maximum airfoil section 
depth is limited to 1.62 inches. Finally, to facilitate manufacture, 
reduce vulnerability to accidental damage, and maintain balance of the 
laminate, the minimum composite skin thickness has been set at five 
plies (0.030 in.). : 


_ Sufficient stiffness of the wing is required to preclude static and 
dynamic aeroelastic instabilities throughout the flight envelope. A 
discussion and analyses of these effects are contained in reference (e). 


Stiffness and strength properties of the unidirectional graphite- 
epoxy material and of the several laminate constructions used in the 
wing, as well as those of the core materials are listed in Table I. 
Since, in general, advanced composites exhibit little static yielding, 
strength calculations for the laminated skins are based solely on ulti- 
mate stresses. For the unmanned aircraft, ultimate loads are 25 per- 
cent over Limit. : 


DESCRIPTION OF WING STRUCTURE 


The principal features of the composite wing construction are shown 
in Figure 2. Laminated graphite-epoxy skins, varying in thickness from 
5 to 30 plies (.030 to .180 in.), are adhesive bonded to the alumimm 
honeycomb core. The laminate construction varies over the planforn, 
The number of plies of each orientation is varied to meet local stiff- 
ness and strength requirements and to maintain balanced construction. 
In the center section and in the outboard two-thirds of the exposed 
span, the core density is 4.5 1b./ft.”. Inboard, where shear loads are 
higher, and particularly in the vicinity of the attachments, higher 
density core materials (6.1, 8.1 and 23.0 1b./ft.3) are used. The 
leading edge is a molded solid section of chopped fiber covered by a 
four-ply, £45° laminated skin. 


To close out the center section and the fuselage interface, and to 
distribute the attachment bolt reaction loads, channel sections of +45° 
construction are placed in the core, Attachment bolt loads are trans- 
ferred to the adjoining structure by means of titanium flanges bonded 
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inside the skins and spool-shaped aluminum fittings inserted into the 
core. Across the aft end of the center section, where the bending mo- 
ment is highest and the attachment bolt load is most critical, the 
graphite-epoxy laminated skins are replaced by titanium plates, step- 
lap bonded to the adjoining laminates, Finally to accommodate elec- 
trical wiring for the wing tip antenna, a fiberglass conduit is 
enclosed within the core between the forward edge of the center section 
: and the wing tip, where mounting holes for the tip antenna pod are 

oe provided. 


FINITE ELEMENT MODELLING AND ANALYSIS 


A diagram of the finite element model assembled for analysis using 
NASTRAN is shown in Figure 3. Triangular and quadrilateral orthotropic 
plate elements are used to represent the skins and core. The stiffness 
(bending and transverse shear) and mass properties of each element are 
derived from the local laminate construction, mmber of plies, core 
properties, and the airfoil thickness, The remaining structural ele- 
ments - leading edge, conduit, channel section ribs and tip antenna pod 
are represented by bar elements. In order to accurately simulate the 
mass distribution of the wing for the vibration analysis, additional 
lumped masses are placed along the leading edge and tip. A listing of 
the NASTRAN bulk data deck is reproduced in Appendix A. 


For the static analysis, the test loads from reference (h) were 
distributed over the finite element model as concentrated forces 
applied to the grid points. To simulate the boundary condition for the 
critical load condition (5g symmetric pull-up), vertical disp Lacements 
were constrained for the grid points located at the fuselage attachment 
bolts, and rotations about the longitudinal axis were constrained along 
the aircraft center line. The results of the static analysis, including 
the deflected shape, deflections at selected points, and reaction loads 
are shown in Figure 4 and Tables II and III. A complete Listing of the 
NASTRAN output data is reproduced in Appendix B. 


To provide vibration frequency and mode data for the modal flutter 

analysis described in reference (e), the real eigenvalue analysis 

option available in NASTRAN was used, Both symmetric and antisymmetric 
es modes were computed by altering the constraints at the fuselage center- 
a line. Similar frequencies, generalized masses and generalized stiffness 
resulted for the two cases. These data are listed in Table IV for the 
first five modes of each case, The mode shapes of the first three 
symmetric modes are shown in Figures 5, 6 and 7, and of the first three 
antisymmetric modes in Figures 8, 9 and 10. A complete listing of the 
NASTRAN output data for the symmetric modes is reproduced in Appendix C. 


NADC-73235-30 


STRESS ANALYSIS 


fo determine the ply stresses in the graphite-epoxy skins, ele- 
ment stresses computed by NASTRAN are first transformed into the local 
laminate coordinates. The stresses in laminate coordinates are shown for 
two streamwise sections in Figures ll and 12. Laminate stiffness prop- 
erties are used to calculate the resulting strains, which are subse- 
quently transformed for each of the ply orientations of the local 
laminate construction. Unidirectional material stiffness coefficients ° 
are then used to compute the ply stresses. To account for the com- 
bined-stress state of the material in determining margins of safety, 
the interaction formula below is applied: 


Oj I 
| x; 7 R; U{t. Mio: a 25% Rae il 
i= 1,2,6 
where Of, , 02, % are the inplane normal and shear stresses 


in the ply and Xj; are the allowable stresses. The margin of safety 
of the element is the lowest of the margins for the individual ply 
orientations. The critical elements of the finite element model and 
the stresses and margins of safety resulting from the static analysis 
are listed in Table V. 


Depending upon the ratio of the skin thickness to the section depth 
(0. € t/y<0.5) the core shear stress varies from 1.0 to 1.5 times the 


average shear stress (force/area), but may be conservatively approximated 
by the relation 


oe < + (1. + tts +) x shear force per unit width 


To determine the core shear stresses and margins of safety, the element 
shear forces (per unit width) from NASTRAN and the skin thickness and 
section depth at the element centroid are used in the fornmla above. 
The resulting stresses are transformed into the ribbon and transverse 
directions of the core, and the strength criterion below is applied: 


we ee 
we et 
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a . e " 0" 


where the subscripto L and W refer to the ribbon and transverse di- 
rections respectively. Shear stresses and margins of safety for the 
critical elements of the model are shown in Table VI. 


The minimum ultimate margins of safety calculated for skins and 
core are 0.14 and 0.32 respectively. Therefore, the wing is con- 
sidered safe for flight. 
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CONCLUSIONS 


1. Based.on the static aerodynamic loads prescribed, adequate margins 
of safety have been provided for the skins and core of the graphite- 
epoxy composite wing designed and fabricated for the BQM-34E target 
vehicle. : 


2. Static and dynamic aeroelastic (flutter) analyses, reviewed in 
reference (e), indicate that the stiffness and mass distribution of 

the composite wing are sufficient to avoid any instabilities throughout 
the flight envelope. 
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TABLE I 


(a) Composite Material and Laminate Properties 


Laminate Construction* Stiffness Coeficients psi x 10° 

L M N of Q11 212 Q52 266 

1 0 0 - 22.12 0.386 1.21 0.40 

1 0 4 45° 9,07 3.99 5.30 3.94 

2 0 4 45° 10.93 3239 4.65 3.70 

4 0 6 45° 11.86 3,09 4.32 3.40 
68 0 2 45° 17.43 1.29 2.34 1.60 

2 0 4 22%° 16.88 1,89 1.70 2.20 


*L = No. of 0° plies; M =™ No. of 90° plies; N * No. of +& plies. 


Unidirectional material allowable ultimate stresses: 


2 


X, = 81.0 ksi X, * 3.6 ksi Xe = 4,05 ksi 


based on average results of specimen tests reduced by 55%: 20% for sta- 
tistical variation; 20% for possible required repair; and 15% for environ- 


mental degradation. 


(b) Honeycomb Core Properties 
Density, lb. /£t.7 Shear Modulus, ksi 


Long. Transv. 
4.5 70 28 
6.1 102 38 
8.1 143 51 
8 


Shear Strength, psi 


Long. Transv. 

350 205 : 
525 305 

740 440 


NADC-73235-30 


TABLE II 


STATIC DEFLECTION AT SELECTED POINTS 


GRID PT. DEFLECTION (IN.) 
i 1 0.01 
mn 6 0.87 
10 5.13 
11 -0.01 
18 1.11 
22 5.41 
35 -0.05 
42 1.64 
46 | 5.98 
59 -0.10 
66 2.15 
70 6.51 
81 -0.03 
| 85 1.91 
ve 90 7.16 
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TABLE IIT 
(a) Attachment Bolt Loads 
Grid Pt. Load (1b.) 

13 - 331.2 
25 - 277.0 
37 106.7 
49 2163.5 
61 3204.7 

TOTAL 4866.7 


(b) Centerline Bending Moment 


Grid Pt. Moment (in.-1b.) 
11 } 4129.4 — 
233 12770.5 
35 | 25466.2 
47 28630.0 
59 : 14849.9 
101 18499.0 
TOTAL 104345 ,0 
10 


MODE 
NO, 


(a) Symmetric 
1 


2 


3 


(b) Antisymmetric 
1 


2 


3 


NADC-73235-30 


TABLE IV 


VIBRATION MODE DATA 


FREQUENCY GENERALIZED MASS 


(HZ) (1b.-in.-sec“*) 
20.2 00476 
69.3 00281 
85.6 «00250 | 
135.1 00344 
143.4 . 00160 
20.6 - 00464 
71.5 - 00266 
85.9 00241 
138.6 . 00381 
144.7 . 00206 


ll 


GENERALIZED STIFFNESS 
(lb.-in.) 


76.5 
532.5 
724.7 

2479.2 


1298.8 


77.4 
535.7 
702.9 

2892.3 


1703.3 
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TABLE V 


LAMINATE STRESS DATA AT LIMIT LOAD 


Laminate 

Stresses, 
psi 

Ox Oy 

22228 2290 
22614 2567 
22442 2098 
20732 1636 
20558 1264 
21406 3254 
23369 2933 
23681 2951 
23793 2669 
23329 2571 
22371 2105 
19743 2503 
21514 3554 


Crit. 
Ply 
Orient. 


12 


Stresses in 
Critical Ply, 


on 
46826 
47277 
47604 
44440 
44615 
49415 
48424 
49103 
49789 
48888 
47433 
46473 
49154 


psi 


O; 
-379 
-311 
~450 
-511 


~619 


-267 


-233 
~243 
~337 
-345 
-444 
-397 


-183 


: NADC-73235-30 
a 
| 
z 
TABLE VI 
ee CORE SHEAR STRESS DATA AT LIMIT LOAD 
Core Density Shear Stresses, psi Vit. M.S. 
El. No. 1b./ft.3 Ribbon Dir. Transverse 

51 4.5 83. 47. 1.40 

52 4.5 102. 28. 1.38 

68 6.1 99, 157. 0.49 

69 8.1 284. 169. 0.32 

70 6.1 164, 46. 1.20 

73 4.5 110, 33% 1.16 

74 4.5 113. 18. 1.33 

75 4.5 111. 2l. 1.33 

92 6.1 167. 27. 1.35 

700 6.1 173. 94, 0.65 

900 8.1 206. 116. 0.99 
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FIGURE 1. 3BQM-34E METAL WING PLANFORM 
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FIGURE 2. COMPOSITE WING CONSTRUCTION 
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STRESS IN LAMINATE COORDINATES, KSI 
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FIGURE 11. LAMINATE STRESSES AT Yy = 24 
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FIGURE 12. LAMINATE STRESSES AT Yy = 29 
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APPENDIX A 


NASTRAN BULK DATA DECK 
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